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Abstract 



An approximation algorithm is proposed to transform truncated QCD 
(or QED) series for observables. The approximation is a modification of 
the Baker-Gammel approximants, and is independent of the renormalization 
scale (RScl) fi - the coupling parameter a(fi) in the series and in the result- 
ing approximants can evolve according to the perturbative renormalization 
group equation (RGE) to any chosen loop order. The proposed algorithm is 
a natural generalization of the recently proposed method of diagonal Pade 
approximants, the latter making the result RScPinvariant in large-/?o approx- 
imation for a(fi). The algorithm described below can extract large amount of 
information from a calculated available truncated perturbative series for an 
observable, by implicitly resumming large classes of diagrams. 
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I. INTRODUCTION 



Pade approximants (PA's) can be regarded as an improvement to truncated perturbative 
series (TPS's), since they have the same formal accuracy as the latter when expanded in 
powers of the perturbative coupling parameter. In addition, PA's act as a kind of analytical 
continuation to TPS's, and thus contain explicitly some additional information not explicitly 
(but implicitly) contained in the TPS's. PA's have been shown to have interesting applica- 
tions in statistical physics and quantum field theory (QFT) [|I|] and in QCD ||. In the latter 
works P|, comparisons of the PA method with other methods which look for optimization of 
the TPS result via a judicial choice of the renormalization scale (RScl) and scheme (RSch) 
were performed and showed good numerical agreement. These other methods include the 
principle of minimal sensitivity 0, the BLM approach 0] and its extensions ||- [[/J, and 
the effective charge approach ||. A novel RSch-invariant method, based partly on the effec- 
tive charge approach, has recently been developed |§, and its precise relations with the PA 
methods, as well as with the method presented here, remain to be investigated. Another new 
approach |K| reduces the RScl- and RSch- dependence by a method of analytic continuation 



from the Euclidean to the time-like region of a kinetic variable. In this context, we mention 
also a review of the role of power expansions in QFT [|ll|] . 



Recently, Gardi (T^j noted that the diagonal Pade approximants (dPA's) to truncated 
series for observables are invariant under the change of the renormalization scale (RScl) \i 
in the large-/? approximation, i.e., when the QCD coupling parameter a(/z) is assumed to 
evolve according to the one- loop renormalization group equation (RGE). His observation 
was based on the invariance of dPA's under the homographic transformations of the dPA 
argument (Ref. Part I): z \— > az/(l + bz). Since the (full knowledge of) observables 



must be RScl-independent, Gardi's observation strongly suggested that the dPA method 
correctly sums up certain classes of multi-loop Feynman diagrams. Since the obtained dPA 
result is invariant in the large-/3o approximation, it was conjectured that the summed-up 
diagrams are one-gluon exchange diagrams where the gluon propagator contains bubble- 



type of radiative corrections. Later, it was explicitly shown [14] that the latter conjecture in 



terms of Feynman diagrams was correct. The authors of |nj stressed that an extension of the 
approximation beyond the large-/?o approximation, i.e., beyond the dPA's, remains a major 
outstanding question in the strive to extract large amount of possible physical information 
from a limited number of available (calculated) terms in the formal QCD perturbation series 
of an observable. 

In this paper, we present an algorithm which goes beyond the large-/5o approximation. 
It it based on a modified version of the diagonal Baker-Gammel approximants (dBGA's)[| 
and results in an expression which has the following properties: 

1. it reproduces the given truncated series up to any even order; 

2. it is invariant under the change of the RScl fi, where the QCD coupling parameter 
a(/i) appearing in the expression can be taken to evolve under the RGE to any chosen 



1 For the conventional Baker-Gammel approximants, see for example Part II of Ref. [[0|. See also 
later discussion. 
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loop order. 

The presented algorithm allows one to extract a large amount of physical information - 
in the sense of the mentioned two points - from a given QCD TPS for any specific QCD 
observable if the latter is available up to (and including) an even order. However, it should 
be mentioned that the approximants of the approach of the present paper, being relatively 
closely related to the usual PA's, may well not be able to discern in QCD nonperturbative 
behavior arising from (UV and IR) renormalons - see discussions in Refs. ||, fl2fl , [14| for 
the case of PA's. The latter behavior is expected to manifest itselfs via a factorial growth of 
coefficients in TPS's. The presented method may prove useful also in areas of physics other 
than QCD or QED. 



II. GENERAL METHOD 



By scaling and raising to a power, it is possible to redefine a generic QCD (or QED) 
observable S in such a way that its (formal) perturbative series acquires the following di- 
mensionless form 



S = Zl f^(z 1 )=z 1 



l + r?>* 



(1) 



where z\ = a(Q\)/n, Q\ is a chosen renormalization scale (RScl), and superscript in coef- 
ficients rW implies that they also depend on the RScl Q\: = r n {Q\). Since the RScl is 
chosen arbitrarily, the full result S is RScl-independent.Q However, as a rule, only a very lim- 



ited number of coefficients A 1 ' are available (calculated), i.e., only a truncated perturbation 



series 



1 +r i { ) z l + r^z\ + ■ 



(2) 



is available, and it is therefore explicitly RScl-dependent as indicated by the superscripts: 
= S n (Ql). Hence, there appears naturally the question of how to construct, on the 
basis of (fj), an expression which satisfies the two previously mentioned points. We will 
now present an algorithm for constructing just such approximations to S'W's (and therefore 
to S), and subsequently show with two theorems and their proofs that the constructed 
approximations satisfy the two mentioned points. 

The evolving QCD (or QED) coupling parameter a satisfies the following RGE: 



dx 



j=0 



i+2 



X 



a(p 2 )/n, dt = <ilnp z 



(3) 



where j3j are constants (RScl-independent), and (3q and j3i are even renormalizat ion-scheme- 
independent (RSch-independent). We can choose to take into account any number of terms 
on the right of this RGE. Now we define the quantity 



2 However, S does depend on a certain energy scale Q typical for the process involved. In scattering 
processes, Q is of the order of the center-of-mass energy. 
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ot(p 2 ) 

fc(Zl,«l) = ; , 2 , 



Zl = a(Ql)/n, Mi = ln(p 2 /Q?) 



(4) 



Its formal Taylor expansion in powers of Ui=ln(p 2 /Qf) is 

00 • 1 d j 

k{z x ,u x ) = 1 + yV^Zi) , where: fcj-(zi) = -77—7 



fc (^i 5 w)L= 



On the basis of RGE we have 

^•(^o = i-iypizi + o{z{ +i ) , fct,^) = 1 . 



(5) 



(6) 



Now we simply rearrange the formal series ([!]) for S/zi, which is in powers of zi = a(Q\)/7r, 
into a related series in kj(zi) 



S = zi 



(7) 



The superscript in the coefficients now again means that they are functions of the RScl 

h 



Q\: fj=fj(Qf). We now define the corresponding formal series in powers of (—Zi) 



^ (1) (^i) 

and construct for z\J-( x \z\) the diagonal Pade approximants (dPA's) 



(8) 



Zl [M -l/M} rW { Zl ) = Zl 



M-l 



i+E 



Q"mZ\ 



m=l 



M 



-, -1 



1 + e w 



n=l 



z 1 ^ 1 \z 1 ) = zi[M — l/M] rm ( Zl ) + O (zf M+1 ) . 



(9) 
(10) 



We recall that & ] is a function of only r\ ^ , 



the above dPA, depending only on / 



(i) 



,rj , due to relations (||). Therefore, also 



'i = h 



2M— 1) by fllOD , depends only on the first 



2M-1 coefficients rf } (i = l,. . . ,2M-1) of the original series ([!]). At this point of algorithm, 
we perform decomposition of the dPA (191) into simple fractions 



z 1 [M-l/M] rm (z 1 )=z 1 ^2- 

i=l v 1 



a. 



+ Wi^l 



(11) 



where [— 1/ui] are the M zeros of the denominator polynomial of the dPA @. As long as the 
latter polynomial doesn't have multiple zeros,Q this decomposition is possible and unique. 
Although the sum ( |TTj) is real for real z%, parameters Ui and Sj are in general complex 



We assume that the exceptional situation of multiple zeros of the denominator does not appear. 
See also later discussion following Eq. (|37|). 
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numbers. They depend on the RScl Q\, because they are functions of r[ ,. . . ^m-i which 
themselves are functions of Q\. Now we come to the central part of the algorithm, by 
constructing a modified version of the diagonal Baker-Gammel approximants (dBGA's) to 
the observable S = z 1 f^ of Eq. ([J) 



M 



ziG [ ^ n 1/M \zi) = zi^2acik(zi,Ui) . 



(12) 



4 = 1 



Function k(z, u), which depends on two arguments and is in our case defined via (^) and ([$]), 
can be called the kernel function of the above modified dBGA. As mentioned previously, the 
dBGA ( |i~2"D is uniquely determined by the first (2M— 1) coefficients r$ (i = \ y . . ,2M— 1) of 
the original perturbative series (D), i.e., knowledge of the truncated series SW-i (@) uniquely 
determines fll2]) . Now we show, by proving two theorems, that these modified dBGA's have 
precisely the two wanted properties mentioned earlier on. 

Theorem 1 ("Approximation" theorem): The modified dBGA ([12]) of order (2M — 1) 
approximates the observable S = Zxp x \zi) of Eq. (|I|) up to (and including) 0(zf M ): 



s = z 1 g%; 1/m \z 1 ) + o(z! m + 1 



(13) 



Proof: We use for k(zi,Ui) in the definition ( |T2"D of the modified dBGA the formal Taylor 
expansion (Bf) in powers of Ui 



M 



ZiG [ ™ n 1/M] (z 1 ) = zi^2&i ufk^zx) = z\ ^ k m (zi) 



i=l m=0 



m=0 



M 



i=l 



(14) 



Here we implicitly assume convergence of the above series. We do the same kind of formal 
expansion for the corresponding dPA (p|)-(|TID 



M 



Zl [M -l/M]f W ( Zl ) 



r m 



i=l m=0 



m=0 



.1=1 



(15) 



According to (0), this expression reproduces ("approximates") the series (||) up to, and 
including, the term Zi/2m_i( — Z\) 2M ~ l ^ z\ M . Therefore, 



M 



(i) 



(m = 0,l,...,2M-l) . 



(16) 



i=l 



Using this in fll4]) and comparing with the full formal series (|7|) for the observable S 
zif^(zi), we obtain 



2M-1 



z 1 G [ ™; 1/M \z 1 ) = z 1 £ f^k m (z 1 ) + 0(z 1 k 2M (z l )) 



O y r [M-l/M] 



m=0 

oo 

*l)=*l E f£ ) hr^l) = 0{z 1 k 2M {z 1 )) = 0(zl M+1 ) 
m=2M 



(17) 
(18) 



where at the end we used estimates @ [kj(zi) ~z%\. Relation ( ]T~8|) proves the theorem. We 
stress, however, that the formal series involved in the proof were implicitly assumed to be 
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convergent. The question of when the proof survives once we abandon the assumption of 
convergence is left open. 

Theorem 2 (Invariance under argument transformation): The modified diagonal Baker- 
Gammel approximant (dBGA) ziG^l^ IM \z\) , as defined in Eq. ([12]), for the observable 
S = zif( l \zi) of Eq. (0), is invariant under the following transformations of the argument: 

Z\ i — ► z 2 (zi]U2i) such that: Z\k{z\,u\) = z 2 k(z 2 ,U\ — U21) , (19) 

where U21 is any (arbitrary) fixed complex number, u\ is arbitrary complex number, 
&i{z\\ W21) is independent of u\, and the following relations are assumed for the fc-function 
(and its Taylor coefficients) appearing in the dBGA: 



k(z,u) ~ z°(— 1) , kj(z) ~ z-* (=>■ ~ Z2) 



(20) 



iVoto to Theorem 2: The seemingly artificial and complicated set of conditions for ar- 
gument transformation is motivated by the following QCD (or QED) interpretation of the 
above parameters: 



a{Q\) 



71 



Ml 



P 



a(p 2 



(21) 



We see then from here and from (|B|) that conditions fl2"U| ) for Theorem 2 are fulfilled in the 
QCD (or QED) case and that the argument transformation in this case means simply the 
change of the renormalization scale (RScl) from Qi to Q2 



Zi = 



7T 



2 2 



a(Ql) 



7T 



[Ul - U 2 l, 



U2j 



k(z 2 ,u 2 ) 



a(p 2 



I) ' 



(22) 



thus implying RScl-invariance of the dBGA ([12]). 

Proof of Theorem 2: Let z\ and z 2 be related by transformation flTD|) . We have two 
different formal series of S in kj(z±) and kj(z 2 ) 



S (= Si/a 



«1 



[In QCD or QED: /• =fj(Q\) and fj Z> =fj(Q2)] The corresponding two formal series in 
simple powers of z\ and z 2 , respectively, are defined in the described algorithm by M) 



c(2) 



22 



(=^/ (2) (^i)) • (23) 



Zl rW(z 1 ) = z 1 



22 



22 



2 2 



1 + E/; 



( 2 ) ( 



•22 



(24) 



Following the algorithm, we construct the usual dPA's for Z\!F^\zi) and z 2 J r ^(z 2 ) in the 
form of decomposition into simple fractions ( |Tl~l) 



M 



2l 



[M-1/M]^i, = S 77 

i=l v 1 



+ u i z lj 



z 2 [M-l/M}^ ) (z 2 ) = z 2 J2 



1 (1 + M^ 2 ) 



(25) 
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again assuming that the exceptional case of multiple poles doesn't appear. The modified 
dBGA's are then constructed in both cases according to the algorithm, by Eq. (|12| ) 



M 



ziG [ m ' J (2i) = zx^2 oi i k(z 1 ,u i ) , 



Z 2 G% } 1/M] (Z 2 ) = Z 2 J2 a M Z 2; 



M 



i=l 



'/(2) 



(26) 



We will now show that ai = a,i and Ui = Ui — u%\. Define expression 



z 2 Q 



[M-l/M] 

/(2) I 



M 



z 2 ) = z 2 &ik(z 2 , Ui - u 21 ) . 



Transformation (p~9|) Zi i— > 22 and definitions (|26|) and (p^ ) imply 



[M-l/M] 

/( 2) I 



^2) — ziGtn 



[M-l/M], 



/(i) 



*1 



Therefore, by Theorem 1 [Eq. ([13])] we have 



S - z 2 Q 



[M-l/M] 

/(2) 



(^) = 5 - ^GjS,- 1 ^!) ~ 



2M+1 2M+1 



Zo 



z 2 k 2M (z 2 ) 



(27) 



(28) 



(29) 



where at the end we used conditions (EDI). We now make the formal Taylor expansion of 
( p7|) in powers of (ttj — w 2 i) 



[M-l/M] 

/(2) 1 



M 00 



Z2) = z 2 Y Oii ^ ( u i - u 2i) m k m (z 2 ) = z 2 Y k m (z 2 ) 



i=l m=0 



m=0 



' M 

E 

i=l 



Qi [Ui - U 2 1 



(30) 



Comparing ( ^9|) and (|30|) , and using (|3|) for the case of argument z 2 , we obtain relations 
analogous to those for the case of argument z\ fll6|) 



M 



(«i - w 2 i) 

8=1 



= (m = 0,l,...,2M-l) 

Therefore, the dPA for z 2 -^ (2) (^) of (H is 



(31) 



M 



z 2 [M-1/M] :f(2) (z 2 )=z 2 J2t [ 



(ui - U21) z 2 ] 



Z2 E (-^2) 



m=0 



Ai 



J2a>i - u 21 ) T 



(32) 



since this dPA expression reproduces the power terms of z 2 T^ (z 2 ) up to (and including) 
~^| M , as can be seen from explicit expansions fl24]) , Cj32|) and relations fl3~l|) 



M 



^ (2) (^)-^Ett- 

1=1 I 1 



+ («i -tt2l)^2 



O (z 2 2M+1 ) 



[M-l/M], 

/(2) 



Z2 



(33) 



to the 



Therefore, this shows that (p7j) is in fact the modified dBGA 

observable S 1 for the case of argument z 2 , according to the described algorithm of constructing 
the modified dBGA's (0) 



M 



Z 2 ) 



Z 2 ^ ®M Z 2, Ui - U 21 ) = Z 2 G 



i=l 



[M-l/M] 

/(2) 



*2 • 



(34) 
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This, together with equality fl28|), shows that the modified dBGA of order (2M — 1), as 
defined by the described algorithm leading to (|l"2"D , is really invariant under transformations 
of argument fllTf). Thus the proof of the theorem is completed. 

Again, it should be emphasized that it is left open under which circumstances the proof 
remains valid once we abandon the assumption of convergence of the series involved. 

In the described approach, the special QCD (or QED) case of one-loop evolution of a{p 2 ) 
(large-/9o approximation) means: k(z\, u-y) — l/(l+floZiUi), kj(z) = /%(— z±Y . Therefore, in 
the one-loop case, expansion ([7|) for S/z\ and (§) for ^^(zi) are identical if Z\ in the latter 
series is replaced by /3 Zi (rescaling). The dBGA (|T^j) is in this case reduced to the usual 
dPA (OJ), with Z\ in the denominators replaced by $qZ\. 

The quantity S/z x = of (0) can lead, at least in some cases, to nonmodified 

dBGA's, i.e., those for which the series (0) has the special property of being a Hamburger 
Part II). This series is the formal Taylor expansion of a Hamburger 



series (cf. Ref. 
function 
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1 + U\Z\ 



where <$>{u\) is increasing with increasing u±, and the coefficients in 



moments 



(i) 



u{d4>{ui) (j 



0,1,2,...) , 



(/■) 



(i) 



(35) 

are the finite 
(36) 



In this special case, observable S of ([I]) has the integral representation 

roc 

S = z 1 p\z 1 )=l zik{z-]_,u\)d^{u\) . 



11) of the dPA's are real 



(37) 



In such (special) cases, the parameters in decomposition 
are positive) |TJJ, and the resulting dBGA's (|T2]) are therefore manifestly real numbers. 
Decomposition (|TT| ) is in this case always possible, i.e., the dPA (Qj has no multiple zeros. 
The case of ([37]), reinterpreted in terms of QCD parameters fl2T|) , means that in the integral 
for S (over du\ = dp 2 /p 2 ) there is a positive function p{u\) = 4>'{u\). The latter can be 
interpreted as momentum probability distribution (see also || and Ul4| ) in diagrams 
with exchange of one effective virtual gluon, where the gluonic propagator contains radiative 
corrections. In fact, if the (approximate) S contains only contributions which can be reduced 
to such classes of one-gluon-exchange diagrams, the authors of |TJJ conjectured that function 
p(ui) is then positive definite. Their conjecture would therefore suggest the following: in the 
cases where the truncated series of @ contains only contributions of diagrams with at 
most one gluon exchange in each one-particle-irreducible part, the method of the modified 
dBGA's presented here would result in a Hamburger series for of (§) and hence in 
nonmodified dBGA's, i.e., with <5j and Ui in (O) being real (a* > 0). Having a truncated 



We note that authors of |13| discuss a somewhat more constrained version of nonmodified 
(d)BGA's, i.e., the one for which ^^(^i) is a Stieltjes series, corresponding to expression (35) 
with lower integration bound being zero. 
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series S^' of (§), with n = 2M — 1 odd, it is straightforward to check whether such a case 
sets in (cf. Part I of |13|| ): determinants of a limited set of matrices A[0],. . . ,A[(n— 1)/2] have 
to be positive, where A[m] is a (m + 1) x (m + 1) matrix whose elements are -Ay[m] =/^-_ 2 
(z,j = l,. . . ,m+ 1). For example, for n—1 (S[ ) this case always sets in, and for n = 3 (S^) 
it sets in whenQ / 2 (1) -(A (1) ) 2 >0. 

We decided to use for the result ( |T2| ) the terminology "modified diagonal Baker-Gammel 
approximant" primarily due to the previously mentioned indirect connection of our result 
with the nonmodified BGA's discussed in ||13|| . The latter reference (part II, Chapter 1.2) 
briefly discusses BGA's, for some specific kernels k(z,u), and from the mathematical and 
numerical point of view. Theorem 1 of the present paper is indirectly related to the Conver- 
gence Theorem 1.2.1 of part II of Ref. [PJ for BGA's, while Theorem 2 of the present paper 



(RScl-invariance) has no analog in [13[ . We stress that the presentation of the present paper 
is self-contained in the sense that the reader is not required to know anything about BGA's, 
but should be reasonably familiar with the usual (d)PA's. In the previous paragraph (which 
does not represent an essential part of the paper), some familiarity with the Hamburger 
series was assumed. 

What to do when parameters oii and Ui in the modified dBGA ( |T2] ) are not simultaneously 
real? In that case, modified dBGA could in principle be complex. However, Theorems 1 and 
2 are valid not just for the entire modified dBGA's, but also for their real (and imaginary) 
parts. Since the observable S is real, we then just take the real part of expression (pT2|). In 
fact, since z x and 5* are real, Theorem 1 [Eq. (|13|) 1 shows that the imaginary part of the 
modified dBGA (^i) must be ~zj M+1 or even less. 

III. EXPLICIT EXAMPLES 

It is instructive to obtain explicit formulas for constructing the dBGA's zxG^ 1 (21) 

of (|12|) once we know the truncated perturbative QCD (or QED) series S2M-1 of Eq. @ for 
the practically possible cases of M=l or M = 2. This will be done below. 

The case M = l in QCD (or QED) is the case investigated already by Brodsky, Lepage 
and Mackenzie (BLM) ||. Formalism presented here, under the mentioned QCD (or QED) 
identifications ( pi] ) and (|3|) for the parameters, gives 

= zjWfa) 1) = zi[l + r? Zl ] =► z 1 G$ ] (z 1 ) = z x k{z u -r?/(3 G ) . (38) 

We note that the above approximation, according to notation ([21]) , is in fact the value of 
the coupling parameter at a scale Q 2 

z 1 G%]\z 1 )=a(Q 2 )/n , Q 2 = Q\ exp(-r 1 1) //3 ) • (39) 

We note that a{p 2 ) evolves according to RGE (|3|) where the number of retained terms (loops) 
can be arbitrarily chosen if known. The result z\Gj${z\) and the scale Q 2 are independent 



5 See also explicit discussion of n = l (M = l) and n = 3 (M = 2) cases later on. 
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of the choice of the RScl Q\ (up to the chosen loop-order) according to Theorem 2, while 
S?[= S 1 (Ql)} and S[ 2) [= S x {Ql)\ in general differ for RScl's Q\ ^ Q\. S{ 2) -S{ 1] ~ z\. In 
fact, it is straightforward to check directly that Z\G^[^{zi) = z 2 G^^(z 2 ) } i.e., that Q 2 of 

(j39|) is RScl-independent, because rj 2 -* —r^ = ri(Q 2 ) — ri(Qf) = /3 ln(Q|/Qf ) / as can 
be checked by using @. This result can be formulated also in a more intuitive manner: 
when choosing the RScl equal to Q 2 of © = z(Q 2 ) i-> z(Q 2 ), rf } = n(Q?) i-> ^(Q 2 ), 
= 5i(Q 2 ) ^ 5i(Q 2 )], the next-to-leading (NLO) term [~ z 2 (Q 2 )] in the perturbative 
expansion of S is zero: S = z(Q 2 )[l+0(z 2 (Q 2 ))]. This follows from fl3"9| ) (i.e., implicitly from 
RScl-invariance Theorem 2 for M = l case) and from Theorem 1 [Eq. fllBDl, under inclusion 
of conditions © [zi ~,z(Q 2 )]. 

The case M = 2 in QCD (or QED) is algebraically more involved, but can be worked out 
in a straightforward way using the presented formalism. Expansion functions kj(zi) (El) are 
obtained from RGE @ 

fci(^i) = -pWx - A*? - ft*? - • • • , (40) 
= + (5/2)PoPi4 + • • • > = - • • • • (41) 

Inverting these relations gives 

z\ = —ksfa) + O(h) , z\ = +-U 2 (^i) + !§to(*i) + O(fci) , (42) 
Po Po z Po 

*i = -^h(zi) - §k(*i) - f||| - ^ fcs(*i) + • (43) 

Po Po \ z Po Po/ 

Inserting this into the truncated series S3 of Eq. (fj) (presumed available), results in the 
rearranged truncated series for S3 of the form ([?[), with coefficients 

/l = —a- , H- -Tp r l + o2*2 , (44) 

Po Po Po 

Ai) _ (_50i , M (i) , (l) _ J_ (i) / 45 n 

Having these coefficients, we construct the truncated series for z\T^' (§), and the dPA 
zi[l/2]- F (i)(zi) (|H]|) in its decomposed form (]TT|). The resulting expressions for parameters 
Ui and a,- are 



«2,1 

where: det 



(/a - /i/ 2 ) ± Vdit 2(/ 2 - A 2 ) ~ , (46) 



/ 3 + /i(2/ 2 -3/ 2 )j +4(/ 2 -/ 2 ) 3 , (47) 
5i = (u 2 - fx)/ (u 2 - u x ) , a 2 = 1 - 5i . (48) 

The plus sign in (^) corresponds to w 2 . For simplicity of notation, we omitted the super- 
scripts in the coefficients = /«(Q?)- Note that expressions (f44|) -([45|) should be inserted 
into ([l6|)-(|48|) in order to obtain these parameters explicitly in terms of the original co- 
efficients r t - = Ti{Q\) (i = 1,2,3). Inserting now these parameters into expression fll2"|) 
(M = 2) gives the sought for RScl-invariant approximation to . Several subcases should 
be distinguished (for M — 2): 
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1. When (/ 2 — ff) >0, then: u^a* are real (z = l,2), Ux^ui and0<aj<l. 



2. When (/ 2 -A 2 )<0 and |/ 3 + /i(2/ 1 2 -3/ 2 )| > 2yJ(ft-h)*, then: ^ are real (i = l,2) 
and U\^U2- 



3. When C/2-A 2 ) <0 and |/ 3 + /i(2/ 1 2 -3/ 2 )| < 2^/(/ 1 2 -/ 2 ) 3 , then: ^ are complex, 
generally complex (i = l,2),Q and ux^u^- 

4. When (/ 2 -A 2 ) <0 and |/ 3 + /i(2/ 1 2 -3/ 2 )| = 2^(/?-/ 2 )3 [or when (/ 2 -/ 2 ) =0 and 

then: the system of equations for ttj and Sj is not solvable, i.e., form (]TTD is 
not valid, the dPA @ has a multiple (double) pole. 

5. When (/ 2 — /i)=0 and fz — fli then: «i,aj are real (i = l,2) and Ux—U2 — fx- 

We should note that in QCD, presently available results of perturbative calculations 
include for various observables S the coefficients and r 2 of (|l|), but not yet r 3 . There- 
fore, at this stage, the algorithm described here still cannot be applied for M = 2 in the case 
of QCD observables (r^f-i = r ^ are not available yet) . This contrasts with QED where 
perturbative coefficients r 3 have been obtained for several QED observables. 

One may raise the question of how to construct, e.g. in QCD, an explicit approximate 
expression for the function k(zx,ux) appearing in the modified dBGA's (|12|) , once we go 
beyond the large-/?o approximation. In the case of QCD (or QED), k(zx,Ux) is defined in 
(pl|). For example, the QCD coupling parameter zx = a s {Q\) /tx & t two-loop level and written 
as an expansion in inverse powers of InQ 2 , is given in 



Zx = 



PoHQl/A 2 



1 - 



(3x In (in ■ 



ln^i 
111 A 2 



(49) 



where the neglected terms are of order (3\ ln 2 [ln(Q 2 /A 2 )]/ ln 3 (Q 2 /A 2 ), and A is the QCD 
scale which depends (like /3j's) on the number of effective quark flavors n/, e.g., A(n/ = 5) ~ 
200 MeV, A(nf = 4) rs 280 MeV. It should be stressed that A is renormalization-scheme- 
dependent (RSch-dependent), hence change of RSch is equivalent to change of RScl when 
effects of /3j's (j>2) are neglected. In order to find k(zx,Ux) [note: U\ = \ia(p 2 /Qf)}, we can 
first write f j49|) for the scale p 2 



Zxk(zx,Ux) = 



1 



7T 



/5 [ln(g 2 /A 2 )+ Ml ] 



1 - 



Pi ln 


In 


(»)+«.] 


PI 


ln( 


f ) + "K 





(50) 



b It can be checked that also in this case the modified dBGA (12) is real, because (|46D-(|48|) imply: 
U2 = {ux)* and a 2 = (ai)*. 

7 Note that convention (||) for QCD coefficients here is different from that in Particle Data Book 
[DJ: f3? DB = 4/3 = ll-2n//3; /3f DB = 8/3i = 51-19n//3; /3^ DB = 128/3 2 = 2857-5033n//9+325n^/27. 



Here, n/ is number of effective quark flavors with mass less than the considered scale Qx- The 
form for /3 2 is in MS scheme; /3q and fix are RSch-independent. 
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To obtain an explicit form for k(zi,ui), we express m(Qi/A 2 ) in fl5"0|) by z\ via ([49]). Inverting 
(H^) can be done numerically, for example by iteration 



In 



01 

A 2 



(n+l) 



1 - 



i Q?l (n) 



In ^J- 





m — r 


(0) 


1 




A 2 







(51) 



If we choose to stop already after the first iteration step, we get 



ln 



Qi 

A 2 



(i> 



^ + | ln(te) 



(52) 



k(z u ui) 



as{p 2 ) 



[1 + p z lUl + {Pxlfa)z x Hf3 Zi)\ 



1+ 



[ln (P oZl ) - ln (1 + + (P 1 /P )z 1 HP oZl ))] \ 

(1 + fozmi + (p 1 /(3 )z l \n((3 oZl )) J 



(53) 



This would then be a fc(zi,wi) function beyond large-/? approximation that could be used 
in the RScl-invariant expressions for modified dBGA's (0), with parameters Ui and oii in 
the case of M = 2 determined from the original truncated series SP by Eqs. ©- 



In the latter equations, would have to be set equal to zero, because (|49|)-(|53D represent 
solutions of RGE ([3|) with only f3 and fi\ retained. The RScl-invariance of the obtained 
modified dBGA in M = 2 case would still not be absolutely precise because (|49f) and (|53D 
represent only an approximate solution to truncated (/3 2 = /3 3 = . . . = 0) RGE (g). To make 
RScl-invariance precise, we would have to integrate the latter truncated RGE numerically 
from t = \n(p 2 /Qf) =0 to t = Ui (i — 1,2) to obtain k(zi,Ui), provided we have at RScl Q\ a 
reasonably reliable value of a s (Qf) =ttzx. 

In fact, numerical integration of the RGE @ would always work and give us such 
k(zi, Ui)'s in (^) that RScl-invariance would be precise. In this RGE, we should not neglect 
at least those j3j coefficients which appear as nonzero numbers in the expressions for the 
coefficients ff (j = 1,. . . ,2M-1) of the rearranged series (|?p, e.g. for M = 2 case we should 
not neglect $i and /?2 [cf. Eqs. (^) -(^5|)1. On the other hand, we can take into account in 
the RGE (|3p more than these minimally required coefficients /3j, thus making the result flT2"D 
RScl-invariant at such an improved level. We note, however, that only a limited number of 
perturbative coefficients (3j (namely: /3 0y . . ,/3 3 ) are known in QCD (cf. M], in MS scheme) 
and QED (cf. 0, in MS, MOM and in on-shell schemes). Therefore, the RGE (|) will 
always be truncated at some level. 

An alternative way to obtain k(zi,Ui), for known z\ and Ui, would be to solve (a trun- 
cated) RGE (|) by integrating it analytically from t = t 1 = \\i{Q\/Q\) =0 to t—Ui, beforehand 
Taylor-expanding the integrands on the right around the point t — and repeatedly using the 
same RGE @ in evaluating Taylor coefficients. If we ignore threshold effects for simplicity, 
this procedure leads to the following expression: 

fc(*i, in) = 1 + z x {-(3 Q Ui) + - folk) + zli-Pffi + ^Miu 2 - foUi) + . . . . (54) 
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If we take into account only a limited number of fa coefficients, than (|53 ) may represent 
a convergent, or a well behaved asymptotic, series for small enough parameters Wj. We 
conjecture that the latter (generally complex) parameters have reasonably small absolute 
values if our choice of Q\ is judicious, i.e., if Q\ is close to a typical scale of the process 
considered. 



IV. CONCLUSIONS, OPEN QUESTIONS 

We presented a method of constructing approximants to available truncated perturba- 
tive series (TPS's) of QED or QCD observables. The approximants are modified diagonal 
Baker-Gammel approximants (dBGA's), with the kernel being a ratio of the gauge coupling 
parameters at different scales. We showed that these approximants have two favorable prop- 
erties: they reproduce the TPS to the available order when expanded in powers of the gauge 
coupling parameter, and they are invariant under the change of the renormalization scale 
(RScl). The gauge coupling parameter can be taken to evolve with the scale at any chosen 
loop-order. 

There are several questions that the presented work raises. For example, it remains un- 
clear whether the dPA's could have multiple poles in some exceptional physical cases, 
and if so, what would have to be changed in the presented formalism in such a case. Fur- 
thermore, the formalism offers RScl-invariant approximants to truncated perturbation series 
of (0) only for odd n = 2M — 1. What would be the best approximant in cases of even 
n = 2M (especially in case n = 2)0? This question remains open even in the large- fa ap- 
proximation when BGA's reduce to the usual PA's. The question of when the series of 
modified dBGA's ZiG^^ l ^ M \zi) (M = 1,2,...), defined by ([12]), is convergent or asymp- 



totic also remains open, especially since our proofs of Theorems 1 and 2 implicitly assumed 
convergence of series involved in the proofs. The latter question appears also when dealing 
with the usual dPA's - however, in such lot of progress in answering this question 



has already been achieved - from formal mathematical [13] and from physical/empirical 
point of view 0. We also note that the dBGA result ([12]) becomes at the level M > 2 
in principle RSch-dependent, since parameters in the dBGA involve also fa (and possibly 
higher fa's) coefficients [cf. ( fPl) -(]48|)1 which are explicitly scheme-dependent. Therefore, it 
would be instructive to test numerically RSch-dependence of dBGA ( |12| ) for M > 2 cases. 
We expect that this RSch-dependence (- dependence on fa) would be very weak, at least 
for such choices of RScl Q\ for which Z\ = a s (Ql)/ir is not large {z\ ~0.09), since work of 
Ref. || (third entry - PRD) indicates that this is the case for some PA's when observable 
5* is the effective charge in the Bjorken sum rule. In any case, one of the next natural steps 
would be to study efficiency of the presented formalism in cases of actual available truncated 
perturbative series for QED observables, to compare it numerically with other methods, and 
to investigate which classes of Feynman diagrams the presented method actually sums up. 



8 The presented algorithm probably cannot be trivially extended to the cases of n = 2M, because 
it relies heavily on the decomposition (11) which is valid only for diagonal PA's. 
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Furthermore, the method, with some modifications, may prove useful also in areas of 
physics others than high energy physics. 
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